Dynamics of Antiferromagnets Driven by Spin Current 
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When a spin-polarized current flows through a ferromagnetic (FM) metal, angular momentum is 
transferred to the magnetization via spin transfer torques. In antiferromagnetic (AFM) materials, 
however, the corresponding problem is unsolved. We derive microscopically the dynamics of an 
AFM system driven by spin current generated by an attached FM polarizer, and find that the spin 
current exerts a driving force on the local staggered order parameter. The mechanism does not rely 
on the conservation of spin angular momentum, nor does it depend on the induced FM moments 
on top the AFM background. Two examples are studied: (i) A domain wall is accelerated to a 
terminal velocity by purely adiabatic effect and no Walker's break-down exists, (ii) Injection of spin 
imbalance significantly splits the spin wave spectrum; a spin current triggers spin wave instability 
above a threshold that depends on the ratio between the spin imbalance and anisotropy gap. 
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Introduction. — Mutual dependence of current and 
magnetization is the central problem of spintronics [1], 
which forms a complementary picture. In a ferromagnetic 
(FM) material where local magnetization varies slowly 
over space and time, conduction electron spins will follow 
the orientation of the background magnetization, known 
as the adiabatic limit. In turn, spin angular momen- 
tum is transferred to the background magnetization via 
spin transfer torques [2-7] as a result of the conservation 
of spin angular momentum. Spin transfer torques pro- 
vide key mechanisms to many intriguing phenomena in 
FM materials such as current-driven domain wall dynam- 
ics [8, 9], spin wave excitations [10-12], etc.. However, in 
antiferromagnetic (AFM) materials, corresponding issues 
are unsolved puzzles hindered by two fundamental diffi- 
culties: (i) staggered AFM order does not respect spin 
conservation with conduction electrons; (ii) adiabatic as- 
sumption is apparently broken by the anti-parallel neigh- 
boring magnetic moments, thus the adiabatic electron 
dynamics in FM materials does not apply. 

On the other hand, many recent experiments [13-15] 
and numerical simulations [16-18] indicate that AFM 
materials exhibit current-induced effects with similar or- 
ders of magnitude, if not stronger than, as those in fer- 
romagnets. Those pioneering investigations ushered the 
field of AFM spintronics [19] and propelled AFM mate- 
rials as promising candidates for real applications. From 
a theoretical point of view, AFM dynamics driven by 
charge current has been studied both phenomenologi- 
cally [20, 21] and microscopically [22, 23]. In the former, 
both adiabatic torque by ac current and non-adiabatic 
torque by dc current are predicted, but an adiabatic ef- 
fect in dc limit is absent; in the latter, adiabatic torque 
can be produced by dc current, but the result includes 
only second order derivatives in space and time. Case be- 
comes rather unclear when turning to spin current, which 
is generated by attaching an AFM metal to a FM polar- 
izer. First of all, it has only been explored phenomeno- 



logically [24] and no microscopic study is yet available. 
But even in the phenomenological model, it is the in- 
duced FM moments on top of the AFM background that 
responds to the spin current, which is a higher order ef- 
fect that drives the AFM staggered order indirectly. Will 
a spin current respond to and drive the staggered order 
directly without the help of induced FM moments? 

We have answered part of this question in a previous 
publication [25] , where adiabatic dynamics of conduction 
electrons is developed in an AFM material with given 
background profile. In this Letter, we solve its converse 
problem — how a spin current exerts back-action on the 
AFM background. In analogy to ferromagnets, electron 
dynamics is adiabatic when the AFM staggered order pa- 
rameter n = rriA — tub {ttia and ttib are neighboring 
magnetic moments) is slowly varying [25]. However, in- 
stead of following the background strictly, electrons are 
subject to internal dynamics between degenerate bands, 
which results in mistracking with the background even 
in the adiabatic limit. The underlying physics is that 
the anti-parallel moments introduce an internal degree 
of freedom on conduction electrons which absorbs dy- 
namics within a unit cell, while dynamics between unit 
cells is governed by the slowly- varying n{r^t) [26]. We 
will follow the same idea but the target here is to derive 
the equation of motion for n{r,t). 

Formalism. — We adopt the Lagrangian approach 
where the system Lagrangian L is 
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with d being the dimensionality. Cn and Cint represent 
Lagrangian densities of the AFM background and its in- 
teraction with conduction electrons, respectively. When 
n{r^t) is slowly- varying, the AFM background is effec- 
tively described by the non-linear a model [27] 
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where c = 2aSJ/h denotes the spin wave velocity, g = 
2\fda^~^ jhS is the couphng coefficient with a being the 
lattice constant and S being the spin of core magnetic 
moments. The last term in Eq. (2) is due to uniaxial 
anisotropy, where n^ includes components of n perpen- 
dicular to the easy axis. 

The interaction term Cint can be obtained by summing 
over contributions of individual electrons. As shown in 
Ref. [25], an electron is restricted to a doubly degen- 
erate band when the exchange coupling is sufficiently 
strong and n{r^t) is slowly varying. The single electron 
Lagrangian Lg is constructed on this degenerate band, 
which yields the electron motion being affected by the 
non-Abelian Berry gauge potential 

A^ = -[-n^ sin Od^if ^ r2^d^0 ^ Ts cos Od^if], (3) 

where r's are Pauli matrices and /i = {t, r} is space-time 
index. 0{r^t) and (p{r^t) are spherical angles specifying 
the local orientation of n(r^t). The parameter ^ repre- 
sents the spatial overlap of Bloch waves of the two sub- 
bands, which is a function of energy and is typically very 
small [28]. Regarding Lg as a functional of n{r^t)^ the 
variational derivative SLe/Sn gives the reaction of the 
electron on the AFM background 
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where ± distinguishes the two sub-bands, and Ve is the 
band velocity of the electron. 

The interaction term Cint is obtained by adding up in- 
dividual electrons Cint = "^Zx I ^^^^e,x{k)fx{k), where 
^^ represents summation over the two degenerate sub- 
bands. Using Eq. (4), we know [28] 
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where (^p is the value of ^ at Fermi energy, ps is the 
equilibrium spin density injected from a FM polarizer 
coupled to the AFM system (e.^., through exchange bias 
on the interface); js represents spin current density that 
can either be injected from an FM material or induced 
by optical means through spin photovoltaic effect [29]. 
Here ps is chosen to be dimensionless — the number dif- 
ference of opposite spin species per unit cell, and the unit 
oi js is fixed to be m/s. As explained in Ref. [28], both 
coherent dynamics within the degenerate sector and in- 
coherent spin-flip scattering due to impurities are highly 
suppressed by the smallness of ^f, thus AFM metals are 
better spin-preservers compared to normal metals. We 
will assume that the system size is within the spin diffu- 
sion length, thus treating ps and the magnitude of js as 
spatially uniform in the following. 
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FIG. 1: (Color online) Schematic view of a setup of AFM 
DW between two pinning ferromagnets at its ends. DW dy- 
namics is described by two collective coordinates, the center 
position Zc and the out-of-plane angle Lp. The DW width W 
is approximately invariant during the motion. 



The last step towards the equation of motion for n{r,t) 
is to account for Gilbert damping by the Rayleigh's dis- 
sipation function R = J d^rlZ = ajd^rn^, and the full 
Lagrangian equation reads 



dC_6^dC_ dC _ dn 

dn dt dh d{\/n) dh ' 
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In view of the constraint n^ = 1, we obtain from Eq. (6) 
the central result of this Letter, 

n X [d^n — (?V^n + cjqTIx] + olu x dtu 

^G{psdt^js-^)n = 0, (7) 

where G = sa^^ ~ ^f) denotes the coupling strength, 
and a — 2ac^/da^~^ /hS is the effective damping coeffi- 
cient. Gompare Eq. (7) with the Landau-Lifshitz-Gilbert 
equation including spin transfer torques in FM materi- 
als [2-7], three distinct features deserve emphasis: (i) It 
is second order in time derivative of n, thus the term 
js ' Vn, though similar to the adiabatic torque in ferro- 
magnets, does not behave like a torque here, it should 
be interpreted as a driving force, (ii) The equilibrium 
spin density ps affects the dynamics explicitly in Eq. (7), 
whereas in FM metals it only exerts an implicit influ- 
ence on the dynamics by slightly renormalizing Gilbert 
damping a and gyromagnetic ratio 7 [4]. (iii) In Fourier 
transform, the ffist two terms and the last two terms 
both give linear dispersions, but the slopes are different 
— spin wave velocity c and drift velocity of the degen- 
erate band Vs = js/ Psi respectively. Their competition 
greatly modifies the spectrum. 

Domain Wall Dynamics. — Due to the absence of dipo- 
lar interaction, formation of an AFM domain wall (DW) 
requires two pinning ferromagnets (along the easy axis) 
at the ends. The pinning originates from exchange bias 
effect on the interface between FM and AFM materi- 
als [30]. Gonsider a DW of 180 degree depicted in Fig. 1, 
such a texture can be achieved by first growing two pin- 
ning FM layers on a homogeneous AFM metal, then ro- 
tate one of them to the opposite direction. Though not 
in exact agreement with theoretical prediction [31], it 



has been realized experimentally in many different con- 
texts [32, 33]. As a compromise between exchange in- 
teraction and anisotropy, the DW assumes a soliton pro- 
file [31]. When the DW is moving, we describe it by the 
Walker's ansatz [34] 
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where the first equation states that n- vectors at different 
positions are kept coplanar and have a common out-of- 
plane angle. The second equation implies that the DW 
remains a soliton shape except that its width W{t) varies 
with time and that the DW moves as a whole with an 
instantaneous center position Zc{t). Eq. (8) enables us 
to compute the total Lagrangian as a functional of three 
parameters Zc, (p^ and W. When DW velocity is much 
smaller than c and its rotation rate is much lower than 
cjo, it can be shown that W{t) is essentially a constant of 
motion [28]. Hence we are left with only two dynamical 
variables Zc and (p, known as the collective coordinates 
of the system. Not bothering with an overall factor, the 
system Lagrangian is effectively written as 



L = 
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The Rayleigh's dissipation function can be calculated in 
a similar way R = a{z'^/W^W(p'^). After some straight- 
forward algebra, we obtain the equations of motion 
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which can be solve analytically. We now scale time t by 
a~^ and the DW center velocity Zc by Vs = js/ps'- define 
V^^ = ZcPs/js and i = at, then eliminate (p in Eq. (10) 
we obtain the following formula 
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where G = PsQ/a. Eq. (11) represents the dynamics of a 
damped harmonic oscillator driven by a constant force. 
With initial condition Vp^(O) = 0, the solution is 



K, 



G^ - Ge-'[G cos Gi + sin Gi] 
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(12) 



V^^ reaches a ter- 



which is plotted in Fig. 2. As t ^ oo, 
minal value 1/(1 + 1/G^). Usually G is quite large in real 
AFM metals, thus the terminal DW velocity Zc{oo) ~ Vg, 
which is the drift velocity of the degenerate band quite 
independent of the spin polarization injected from the 
pinning ferromagnet. 

We make a numerical estimate for IrMn (or FeMn): 

the lattice constant a is 3.6 ~ 3.8 A, the Gilbert damping 

-1. 




FIG. 2: Scaled DW velocity as a function of time. When G is 
very big, ic(oo) ^ Vs, thus the larger the current, the faster 
the DW moves, no Walker's break- down exists. 



spin S* is 2 ~ 3 and c is roughly 10^ m/s. So even for an 
injected spin polarization P ~ 0.05, we have G ~ 10^. 
According to Eq. (12), to drive the DW up to 100 m/s, 
the required current density is of 10^ A/cm^, which is 
smaller than that in FM materials. 

To close the argument on DW, we provide some re- 
marks. Due to vanishing net magnetization, AFM DW 
motion is not affected by demagnetization (which leads 
to intrinsic pinning) like their FM counterparts, which 
could partly explain why an AFM DW is easier to drive. 
Moreover, getting rid of demagnetization also removes 
the Walker's break-down, and the DW velocity is simply 
proportional to the drift velocity when G is large. This 
means as we increase the current density, arbitrary DW 
velocity can be achieved. Of course Joule heating for- 
bids the blowing up of current density, but for a sizable 
DW velocity, the required current density turns out to be 
smaller than that in FM materials, which helps to avoid 
large Joule heating. In addition, our theory is based 
on adiabatic electron dynamics, thus G{Psdt + js ' V)n 
only includes the adiabatic effect of spin current. While 
only non-adiabatic torque determines the terminal veloc- 
ity of a FM DW [4], the AFM DW here is driven to a 
steady motion by purely adiabatic forcing, the (transfer) 
efficiency of which is usually much higher than that of 
non-adiabatic effects. This is also responsible for why an 
AFM DW is more movable. 

Spin Wave Excitations. — Injecting spin imbalance and 
spin current significantly modifies spin wave excitations 
in AFM metals. Denoting the easy axis by e, we adopt 
the spin wave Ansatz n = e + n^e^^^''^~^'^\ where n^ 
is a small deviation {\n±_\ <C 1) perpendicular to e. It 
worths mentioning that the relative motion between M^ 
and M^ within a unit cell (dynamics oi m = (M^ + 
M^) which respects m ■ n = 0), which seems to have 
been ignored, is in fact resolved into the dynamics of n 
described by Eq. (2) [27]. Using the Ansatz above, the 
spin wave spectrum is easily obtained. 
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c'k')±Psg{uj-Vs-k) = 0, (13) 



rate is similar to typical FM metals a ~ 10 s ; the core where + (— ) refers to the case where the direction of 



A (B) sublattice is pinned along the FM polarizer. First 
consider the macrospin model with k ^ and the system 
processes as a whole, the spectrum is 



■-[±p,g±Jip,g) 



■4C.2] 



(14) 



where the two ± are independent. Plotted in Fig. 3, we 
see that the splitting is proportional to the spin imbal- 
ance AljJ = PsG' In typical AFM metals such as IrMn, 
even with a rather small injected polarization P ~ 0.05 
from the polarizer, psQ is roughly 100 GHz, which is of 
the same order as (or even larger than) the anisotropy 
gap ujQ. Such a significant splitting can be easily mea- 
sured by AFM resonance [35]. 

We also study the general case with finite /c, and solve 
for uj{k) in the longitudinal direction. As current density 
is increased, the imaginary part ofuj{k) becomes positive 
after a threshold, by which damping turns into ampli- 
fying. As a result, spin wave at certain frequencies be- 
come unstable, i.e., magnons are emitted by fast moving 
electrons. By setting Im[a;(/c)] = 0, the threshold drift 
velocity of the degenerate band is obtained [28] 
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where k^ = ujq/c. The threshold reaches a minimum at 
k = ko^ which marks the most unstable mode. For this 
particular mode, the wave length in typical AFM met- 
als (such as IrMn) is Aq ~ 10^ nm. Since Aq is much 
larger than the lattice spacing, our assumption at the 
beginning is guaranteed. Besides /c, the ratio psG/^o 
also plays an important role, it vanishes in the absence 
of spin imbalance, by which the threshold goes to infin- 
ity and no spin wave instability exists. But as mentioned 
above, a spin injection with P ~ 0.05 in IrMn is enough 
to bring this ratio to the order of 1, hence from Eq. (15) 
the corresponding threshold current density for k = ko is 
estimated to be 10^ A/cm^. Furthermore, it is remark- 
able that the damping coefficient a does not appear in 
Eq. (15), though the instability is physically due to the 
overcoming of damping. 




FIG. 3: Spliting of spin wave spectrum (at zero k) due to 
injected spin imbalance. Even for P ^ 0.05, psQ /ojq can be 
of order 1, thus the splitting is comparable to ujq. 



It can be shown that Re[co'(/c)] is also zero at the thresh- 
old [28], which means the spin wave instability starts 
by forming a stationary pattern n{r) ^ ex.-p[ikor] at the 
threshold point, and gradually evolves into propagating 
mode afterwards. When inhomogeneous spatial configu- 
ration is developed, exchange energy of the AFM back- 
ground is increased. Therefore, to sustain this instabil- 
ity, kinetic energy of conduction electrons must be trans- 
ferred continuously to the background. This results in a 
sudden rise of differential resistance dV/dl at the low- 
est threshold (for k = /cq), which is detectable with high 
accuracy using today's technology [11, 12]. 

Finally we compare Eq. (7) with results from existing 
literature. Within the adiabatic limit and linear order 
in Vn, Ref. [20, 21] give non-trivial results only for ac 
current, whereas dc current produces null effect. In con- 
trast, dc spin current in Eq. (7) does generate first order 
effect in the adiabatic limit. Ref. [22] studied dc charge 
current, but the result turns out to be second order terms 
d'^n^V'^n^ and dtVn. Ref. [24] considered spin current, 
but it only couples to the induced FM moments on top of 
the AFM background, which drives the staggered order 
indirectly, and the result is again second order. 

We thank E. Tveten, A. Brataas, A. Qaiumzadeh, M. 
Tsoi, A. MacDonald, and K. Everschor for helpful discus- 
sions. This work is supported by DOE-DMSE, NBRPG, 
NSFC, and the Welch Foundation. 



Electronic address: rancheng@utexas.edu 

I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 
76, 323 (2004) and the reference therein. 
L. Berger, Phys. Rev. B 54, 9353 (1996); J. Slonczewki, 
J. Magn. Magn. Mater. 159, LI (1996). 
Y. B. Bazaliy, B. A. Jones, and S. -C. Zhang, Phys. Rev. 
B 57, R3213 (1998). 

S. Zhang and Z. Li, Phys. Rev. Lett. 93, 127204 (2004). 
D. C. Ralph and M. D. Stiles, J. Magn. Magn. Mater. 
320, 1190 (2008). 

C. H. Wong and Y. Tserkovnyak, Phys. Rev. B 80, 
184411 (2009); Y. Tserkovnyak and C. H. Wong, Phys. 
Rev. B 79, 014402 (2009). 
A. Brataas, A. D. Kent, and H. Ohno, Nature Materials 

II, 372 (2012). 

G. S. D. Beach, M. Tsoi, J. L. Erskine, J. Magn. Magn. 
Mater. 320, 1272 (2008). 

Y. Tserkovnyak, A. Brattas, and G. E. W. Bauer, J. 
Magn. Magn. Mater. 320, 1282 (2008). 
Z. Li and S. Zhang, Phys. Rev. Lett. 92, 207203 (2004). 
Y. Ji, C. L. Chien, and M. D. Stiles, Phys. Rev. Lett. 90, 
106601 (2003). 

M. Tsoi et a/., Phys. Rev. Lett. 80, 4281 (1998); M. Tsoi, 
V. Tsoi, J. Bass, A. G. M. Jansen, and P. Wyder, Phys. 
Rev. Lett. 89, 246803 (2002). 
Z. Wei et al, Phys. Rev. Lett. 98, 116603 (2007). 
S. Urazhdin and N. Anthony, Phys. Rev. Lett. 99, 046602 
(2007). 

[15] B. G. Park, et al, Nat. Mater. 10, 347 (2011). 



[1] 
[2] 

[3: 

[4] 
[5] 

[6: 

[7] 
[8] 
[9 

[lo: 
[11 

[12 



[13: 

[14 



[16] R. Wieser, E. Y. Vedmedenko, and R. Wiesendanger, 
Phys. Rev. Lett. 106, 067204 (2011). 

[17] Y. Xu, S. Wang, and K. Xia, Phys. Rev. Lett. 100, 
226602 (2008). 

[18] T. Jungwirth, et al Phys. Rev. B 83, 035321 (2011). 

[19] A. H. MacDonald and M. Tsoi, Phil. Trans. R. Soc. A 
369, 3098 (2011). 

[20] K. M. D. Hals, Y. Tserkovnyak, and A. Brataas, Phys. 
Rev. Lett. 106, 107206 (2011). 

[21] E. G. Tveten, A. Qaiumzadeh, O. A. Tretiakov, and A. 
Brataas, Phys. Rev. Lett. 110, 127208 (2013). 

[22] A. C. Swaving and R. A. Duine, J. Phys.: Cond. Mat. 
24, 024223 (2012); A. C. Swaving and R. A. Duine, Phys. 
Rev. B 83, 054428 (2011). 

[23] P. M. Haney and A. H. MacDonald, Phys. Rev. Lett. 
100, 196801 (2008); A. S. Nunez, R. A. Duine, P. Haney, 
and A. H. MacDonald, Phys. Rev. B 73, 214426 (2006). 

[24] H. V. Gomonay, R. V. Kunitsyn, and V. M. Loktev, Phys. 
Rev. B, 85, 134446 (2012); H. V. Gomonay and V. M. 
Loktev, Phys. Rev. B, 81, 144427 (2010). 

[25] R. Cheng and Q. Niu, Phys. Rev. B 86, 245118 (2012). 

[26] Note that even the dynamics within a unit cell (inter- 
nal dynamics) is itself adiabatic, which defines a geomet- 
ric mapping from the n-orbit to electron spin orbit (see 
Ref. [25]). Also see the supplementary material [28]. 



[27] F. D. M. Haldane, Phys. Rev. Lett. 50, 1153 (1983); 
ibid, 61, 1029 (1988); E. Fradkin, Field Theories of Con- 
densed Matter Systems^ Addison- Wesley publishing com- 
pany, 1991. 

[28] See the supplementary material through Link: will be 
available upon publication. 

[29] S. M. Young, F. Zheng, and A. M. Rappe, Phys. Rev. 
Lett. 110, 057201 (2013). 

[30] F. Nolting et a/.. Nature 405, 767 (2000); J. -V. Kim and 
R. L. Stamps, Phys. Rev. B 71, 094405 (2005); D. Mauri, 
H. C. Siegmann, P. S. Bagus, and E. Kay, J. Appl. Phys. 
62, 3047 (1987). 

[31] N. Papanicolaou, Phys. Rev. B 51, 15062 (1995); ibid. 
55, 12290 (1997). 

[32] F. Y. Yang and C. L. Chien, Phys. Rev. Lett. 85, 2597 
(2000). 

[33] M. Bode, E. Y. Vedmedenko, K. von Bergmann, A. Ku- 
betzka, P. Ferriani, S. Heinze, and R. Wiesendanger, Nat. 
Mater. 5, 477 (2006); M. Bode et a/.. Nature 447, 190 
(2007); P. Sessi, N. P. Guisinger, J. R. Guest, and M. 
Bode, Phys. Rev. Lett. 103, 167201 (2009). 

[34] N. L. Schryer and L. R. Walker, J. Appl. Phys. 45, 5406 
(1974). 

[35] F. Keffer and C. Kittel, Phys. Rev. 85, 329 (1952). 



